www.gusucode.com > matlab编程NSCT分解 图像融合 各个融合指标评价体系 分解源码程序 > matlab编程NSCT分解 图像融合 各个融合指标评价体系 分解源码程序/NSCT/nsct1.m
function coeffs=nsct(im,n) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Image decomposition by nonsubsampled contourlet transform (NSSC). % This is the iterated filter bank that computes the nonsubsampled % contourlet transform. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Parameteters: nlevels = n ; % Decomposition level %pfilter = 'pkva' ; % Pyramidal filter dfilter = 'dmaxflat7'; %'cd' ; % Directional filter pfiltr= '9-7'; % Nonsubsampled Contourlet decomposition %im=imread('zoneplate.png'); coeffs = nsctdec( double(im), nlevels ,dfilter,pfiltr); % disp( nlevels); disp(dfilter);disp(pfiltr); % % % Display the coefficients % %disp('Displaying the contourlet coefficients...') ; % %shownsct( coeffs ) ; % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Nonsubsampled Contourlet transform (NSSC) reconstruction. % % This is the inverse of nsscdec, i.e. % % imrec = nsscrec(coeffs, dfilter, pfilter); % % would reconstruct imrec = im % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % % Reconstruct image % imrec = nsctrec( coeffs, dfilter, pfiltr ) ; % % % % Show the reconstruction image and the original image % %if 0 % figure; % subplot(1,2,1), imagesc( im, [0, 255] ); % title('Original image' ) ; % colormap(gray); % axis image off; % subplot(1,2,2), imagesc( imrec, [0, 255] ); % title('Reconstructed image' ) ; % colormap(gray); % axis image off; % %end % % mse = sum( sum( (imrec - double(im)).^2 ) ); % mse = mse / prod(size(im)); % % disp( sprintf('The mean square error is: %f', mse ) ); % disp(' ');